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A charged particle circling in a uniform magnetic field and kicked by an electric field is considered. 
Under the assumption of small magnetic field, an iterative map is developed. Comparison between 
the (relativistic) non-radiative case and the (relativistic) radiative case shows that in both cases one 
can observe a stochastic web structure, and that both cases are qualitatively similar. 



I. INTRODUCTION 

Zaslavskii et al |IJ] studied the behavior of particles in 
the wave packet of an electric field in the presence of 
a static magnetic field. For a broad wave packet with 
sufficiently uniform spectrum, one may show that the 
problem can be stated in terms of an electrically kicked 
harmonic oscillator. For rational ratios between the fre- 
quency of the kicking field and Larmor frequency associ- 
ated with the magnetic field the phase space of the sys- 
tem is covered by a mesh of finite thickness; inside the 
filaments of the mesh, the dynamics of the particle is 
stochastic and outside (in the cells of stability), the dy- 
namics is regular. This structure is called a stochastic 
web. It was found that this pattern covers the entire 
phase plane, permitting the particle to diffuse arbitrarily 
far into the region of high energies (a process analogous 
to Arnol'd diffusion (|]). 

Since the stochastic web leads to unbounded energies, 
several authors have considered the corresponding rela- 
tivistic problem. Longcope and Sudan Q studied this 
system (in effectively 1^ dimensions) and found that for 
initial conditions close to the origin of the phase space 
there is a stochastic web, which is bounded in energy, of 
a form quite similar, in the neighborhood of the origin, 
to the non-relativistic case treated by Zaslavskii et al. 
Karimabadi and Angelopoulos j| studied the case of an 
obliquely propagating wave, and showed that under cer- 
tain conditions, particles can be accelerated to unlimited 
energy through an Arnol'd diffusion in two dimensions. 
Since an accelerated charged particle radiates, it is im- 
portant to study the radiative corrections to this motion. 
We shall use the Lorentz-Dirac equation to compute this 
effect. 

We compute solutions to this equation for the case of 
the kicked oscillator. At low velocities, the stochastic 
web found by Zaslavskii et al (!]] occurs; the system dif- 
fuses in the stochastic region to unbounded energy, as 
found by Karimabadi and Angelopoulos j|. The veloc- 
ity of the particle is light speed limited by the dynamical 
equations, in particular, by the suppression of the action 
of the electric field at velocities approaching the velocity 



of light §. 

II. MODEL 

In the present study we will consider a charged par- 
ticle moving in a uniform magnetic field, and kicked by 
an electric field. The effect of relativity, as well as the 
radiation of the particle, will be considered. We restrict 
ourselves to the "on mass shell" case^j. 

The fundamental equation that we use to study radia- 
tion is the Lorentz-Dirac equation R , is 

e 1 
m i M = -x v F^ + 7 m (i'' 1 - -^x^x^x"), (1) 
c c 

where 7o = f = 6.26 x 10~ 24 sec. /i and v are the coor- 
dinates t and x, y, and z (or 0, 1, 2, 3), and the derivative 
is with respect to r, which can be regarded as proper time 
in the "on-mass-shell" case. is the antisymmetric 

electromagnetic tensor. The first term of the right hand 
side of Eq. (|l|) is the relativistic Lorentz force, while the 
second term is the radiation-reaction term. Note that 
the small size of the radiation coefficient (70 <C 1) leads 
to a singular equation that requires special mathemati- 
cal treatment, as well as physical restrictions, as will be 
shown in the succeeding sections. 

Following Zaslavskii et al jjj, the magnetic field is cho- 
sen to be a uniform field in the z direction, and the kick- 
ing electric field is a function of x in the x direction, 

B= (0,0, B ) 

OO 

E(x,t) = (f(x) Y, S(t-nT),0,0). (2) 

n— — 00 

Originally, Zaslavskii et al chose a uniform broad band 
electric field wave packet which can be expanded as an 



"The more general case which includes also the possibility 
"off mass shell" motion, is discussed in a separate article in 
this journal [0 . 



f 



infinite sum of (kicking) ^-functions (loq is the frequency 
of the central harmonic of the wave packet, fco is the wave 
number of the central harmonic, and Aw is the frequency 
between the harmonics of the wave packet) 

oo 

E x — E(x, t) = —Eq sin(fcoa; — tu^t — nAuit), (3) 

n— — oo 

which, for ujq = 0, becomes 

oo 

E(x,t) = -EoTsiniknx) ^ 5{t-nT), (4) 



• u • ■ -z , -2 2 12 

x^Xfj, — x + y — c t 



n— — oo 



where T = Eq. (||) is the generalization of Eq. 

(Q), with an arbitrary function f(x) instead of the sine 
function of Eq. (H). 

In order to introduce a map which connects between 
cycles of integration which start before a kick and end be- 
fore the next electric field kick, one has first to integrate 
over the 6 electric kick and then to integrate the equa- 
tions of motion between the kicks, where only a uniform 
magnetic field is present and there is no electric field up 
to the beginning of the next kick. 



III. THE MOTION OF A CHARGED PARTICLE 
IN A UNIFORM MAGNETIC FIELD 

In the case of a uniform magnetic field from Eq. (|^), 
Eqs. (|l|) reduce to three coupled differential equations, 

1 



Ct = .T = Jo('Xo 7)XqR) 

cr 

X — X\ — —VtX2 + 7o('i'i nXlR) 

y = x 2 = flx 1 + 7o(i' 2 - \± 2 R), 



(5) 



where fl — R = x\ + x\ — x\ , and e = — e (charge 

of the electron). 

Using a complex coordinate ||, u = x + iy, Eqs. (||) 
can be written as, 
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t = 7o(t - -~tR) 
ii = iilu + 7o(u yiiR). 



(6) 



We wish here to constrain the radiatively perturbed mo- 
tion to the "mass shell" ; to do this, it is very convenient 
to use hyperbolic coordinates, 



i = cosh q 

x = c sinh q cos < 

y = c sinh q sin (, 



(7) 



since the "on mass shell" restriction (see ](| for further 
discussion), 



(8) 



is then automatically satisfied. The complex coordinate 
u then becomes, 



u = c sinh qe 



i<j> 



(9) 



and the equations of motion Eqs. m) can be written as, 



q = — 7o</> 2 cosh q sinh q + 70? 
4> = Q + 2jo(f>q coth<7 + 7o</>- 



(10) 
(11) 



As pointed out earlier, in the case of a singular equa- 
tion (such as Eqs. (10) and (|lll)), one has to consider 
physical arguments as well as mathematical arguments. 
There are several suggested methods to avoid the "run 
away electron" problem Q. One of the most frequently 
used, especially useful for scattering problems, assumes 
H that the particle loses all its energy after a sufficiently 
large time, and the run away parts of the solution are 
set to zero; the equations of motion can then be written 
as an integral equation. Another method, which permits 
the study of problems with only non-asymptotic states 
|l0| , uses an iterative singular perturbation integration 
that leads to the stable solution. Since the integration 
we must use is over a finite time, it is impossible to use 
an asymptotic condition (e.g. \v\ — > at r — > 00). The 
Sokolov and Ternov approach Q is more suitable for our 
purpose, since it can be implemented for bounded time 
integration, and the mathematical solution is simple. In 
this method, in the first step, the perturbation terms 
in Eq. (|lTJ) are disregarded (as in the iterative scheme 
of ||Io|l) , and then the resulting equation is exactly in- 
tegrated. In the second step, the singular term on the 
right hand side of Eq. ( |To| ) is not considered; using </>(r) 
from Eq. (pd|), the equation including just the first term 
is integrated. The solution is § 



4> = Vlt 

[3(t) = tanhq 



c 



(12) 



where fto = ^ is the actual (normalized) velocity, and 
is the decay time for the energy of the particle. 



T() 



M' 2 



To get an estimate for the radiation during one cycle, 
we will refer to the maximal uniform magnetic field that 
can achieved today in a laboratory, which O(10T) = 
O(10 5 Gauss). If, for example, Bq = 10 5 Gauss, then 
= 1.76 x 10 12J -, and r = 5 x 10~ 2 sec. Thus, it is 
clear from Eq. ( |l2j j that the particle makes |^to ~ 10 10 
cycles before it decays to ~ of its initial velocity. In other 
words, the energy loss during one cycle in very small, and 



Tt is possible to obtain a short-lived magnetic field of 
O(10 6 Gauss) in the laboratory. 
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since in our problem, the time T between the kicking 
is of the order of the period, the energy loss between 
consecutive kicks is very small. At much higher field 
strengths, a qualitatively different behavior may occur, 
for which the energy loss may be very high during a single 
cycle, and the resulting motion may be modified to reflect 
a nonrelativistic behavior, or it may entirely stop before 
the next kick. This behavior may result in a stochastic 
web of a somewhat different type. 

The time between the kicking is measured according 
to the observed time along the motion, At — T . Thus, 
it is necessary to find the corresponding At. It follows 
from Eq. (g) and Eq. (||) that, 



i = cosh q = 



I -Pie 



(13) 



The solution of Eq. (|T|) can be obtained by a elementary 
integration, 



t = To In 



1 + s/T 



-Ci 



(14) 



where a 
gets, 



Poe T o . After some algebraic operations one 



In Fig. [I] we present the behavior of the function 
At ^to ' ^ can ^ e seen f rom Fig- that in any case, 



At < T: 



(16) 



this implies that the time difference according to the 
proper time r is always less then the time difference ac- 
cording to observed time t. This fact is similar to the well 
known relativistic "time-dilation" phenomenon. The in- 
equality Eq. ( |l6| ) can be shown also by considering the 
two extreme cases, 



T » r => At » T + T ln 



l + ^l^Wo 



— B 2 T 2 

2 T 



(17) 
(18) 



For strong magnetic field (i.e. the case of Eq. fll7p) 
the numerator is less (or equal) then 2 and thus the In 
function will give a negative number; in this case the 
inequality Eq. ( |k|) is clearly achieved. For a low mag- 
netic field Eq. (18) behaves like a parabola; the slope 
for T = is \/l — ftp which is less than one, and thus 
the inequality Eq. (119) is satisfied. As explained above, 
in the present study we confine ourselves to the case of 
very small T, since this time difference is constrained be- 
cause of the limited magnetic field that can achieved in 



the laboratory. Thus, just Eq. 
study. 



Fig. 1 



lSj) is applicable in this 



At/t„ 




FIG. 1. The function At( ^ /to) versus T/r . Four typical 
cases are shown (/3o = 0, fio = 0.7, /So = 0.9, /3o = 1). All 
possible curve lines between the graphs of /3o = and f)o = 1. 
The dashed line indicates the asymptotic behavior of the case 

A) = l- 



IV. DERIVATION OF THE MAP 

In the previous section we have shown the approximate 
solution for a charged particle circling in a uniform mag- 
netic field. However, as mentioned above, the kicking of 
the electric field from Eq. ^ is with respect to the ob- 
served time i, and it is therefore convenient to integrate 
the equations of motion (Eqs. ( |l0| ) and ([ll])) with respect 
to t. 



A. A charged particle in a uniform magnetic field - 
integration with respect to the observed time t 

The arguments that were used in order to obtain Eqs. 
(|l2| ) , as well as the implementation of the chain rule on 
Eqs. (|l0|) and (0) by the use of Eq. @) lead to the 
following equations of motion 



dq 1 

— = Slim q 

at to 

dd> n 



dt cosh q 
A simple integration of Eq. (|l9|) yields, 



tanh ■ 



e Cl e~ 



(19) 
(20) 

(21) 
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where C\ < 0. Substitution of Eq. ( fH| ) in Eq. ( p0| ) gives, 



-p = fttanh( Ci), 



and the solution is 



firo lncosh( C\) 



C 2 



(22) 



(23) 



where C\ = cosh _1 (^-), and G% = tan _1 (y /io) + 
f2roln/3o- Returning to the original x, y coordinates (us- 
ing Eqs. (0)), the equations of motion become, 

dx cos(r2r ln(cosh(i/r — C\)) + Ca) 
dt cosh(t/r — Ci) 

sin(rir ln(cosh(t/r — C\)) + C2) 
dt cosh(i/ro — C\) 



(24) 



It is possible to write Eq. (E4h in more convenient way, 
by using the relation 



cosh(- 



t 

'T 



Ci) 



— cosh — 

PO T 



+ J-0T^ sinh A. (25) 



Eqs. (|24| ) then become 



dx 
~dt 

dt 



[ft ) cosa 



(*)o 



ana 



cosh^ + Vl-$sinh-L 

(%)„«»" 
cosh 4- + Vl-^sinh^' 



where 



a = fh- ln(cosh — (1 + Jl- /Stanh— )) 
to v t 



(26) 



(27) 



It is necessary to integrate Eq. (|26j) since the value of 
x is used in the electric field kicking in Eq. (||) . There 
is no analytical solution to Eq. (26), and a numerical 
integration must, in general, be performed. However, 
since we restrict ourselves to T/tq <C 1, it is possible to 
expand 4| and ^ in a Taylor series and then to integrate. 
The expansion of Eq. (|27j) is 



a « n^Tfit + + \ V /I~^ g ( | _ ^ + 0( ^ ; 

» Z To O Tq Z Tq 

(28) 

The first order expansion (according to t/Vo) of Eqs. (|2 
is 



It 



dy 
dt 



1 



1 + 



dx 



dt 



1 



cos ( _t + -n$-) 



7 



TO 



sin(-t + -QPo 3 -) 
7 2 t 



sin(-< + ^/3 2 -) 
at y 7 2 t 



1 



,t 2 



dy 



dt 



,9- 



1 



-±) cos(-t+-fi^-) 



7 



TO 



(29) 



where 7=1/ ^/ 1 — pg . Since 
approximation, 



1 + 

7 to 



1 



7 To 

1 t 

7 To' 



<C 1, one can use the 



(30) 



Using Eq. (p3), the actual velocities, dx/dt and dy/dt, 
from Eq. (|29|) can be integrated and expressed by elemen- 
tary Fresnel functions. However, the effect of radiation 
is due to the term (which multiplies the sine and 

cosine functions), and the |n/3o|- term is not essential 
since the major offset from the unstable fixed point is 

due to the —t term. 

7 

Under the above assumptions the resulting equations 



— sin(-T) 

dt I n 7 



T 

1 I 

7To, 

dx\ St 

— cos(-T) 
dtj 7 

i-L), 

7T0/ 

— sin -T + ( -p cos -T 

dtA S ^ Wo 7 . 



(31) 



(32) 



One therefore obtains 



~ 7 



dy 



dx 



(33) 



The exponential decay from Eqs. ( f24[ ) are replaced by 
linear decay. 



B. Integration over a 5 electric kick 

As pointed out in the previous section, the Lorentz- 
Dirac equation (Eq. (Q)) is a singular equation. The 
electric field which is used in this paper was expanded 
to a sum of 6 functions (Eq. (g)). In that case, Eq. (|l|) 
can not be integrated over the kick by regular treatment, 
and some approximation for the S function should be 
considered instead. In the present paper we assumed that 
the radiation during the kick is negligible; the radiation 
during the kick will be studied elsewhere. In addition, 
since the integration is over an infinitesimal time interval, 
there is no need to consider the constant magnetic field 
during the kick. 

Under the above assumptions, in the neighborhood of 
the kick, the Lorentz-Dirac equation, Eq. (|l]), can be 
written as, 



d ( dx 
dt\dr 



= f(x) £ 5{t-nT) 



(34) 



+ Xq. 



4 



Integration over the (5 function yields, 



C. Analysis 



dx\ 


f dx 








(dy 


dr L " 


\ dr 



(35) 



were the + sign indicates after the kick, and the - sign, 
before the kick. Following Ref. Q, we chose, 



eET 

f(x) = sin(fca:) = K sin(fcx). 



to 

Using the fact that -£ 
relations 



(36) 



dr 



7^ one obtains the following 



7+ = 




§1 < f 



+ 7 2 



dx\ 1 „, , 

-77 +—f(Xn) 
dt ) _ 7+ 

dt J ' 



(37) 



Up to the derivation of Eqs. (p8|)-(|40|) we have passed 
three stages, (a) the velocities were calculated using an 
exact solution and the position x can be calculated by 
numerical integration (Eqs. (|2^)), (b) the velocities were 
approximated by Eqs. ( p9| ) and the position x by the 
elementary Fresnel function, and (c) the velocities were 
approximated by Eqs. (|3^)-(39) and the position x (Eq. 
(plO|)) was derived using elementary integration. Among 
all the combinations of the solutions for the velocities 
and the position x we choose the combination of the ex- 
act solution of velocities (case (a), Eqs. (p6|)) and the 
exact solution for the position x (stage (c), Eq. (p0[)). 
The value of x was selected from the third stage of our 
derivation since it enters the equation just in the kicking 
term, and there it affects only slightly the phase. Thus, 
one does not expect that this fact fact can change the 
typical behavior of the particle. The above analysis is 
summarized in Fig. |2[ 

Fig. 2 



Returning to the initial charge e — — cq and thus Q — > 
— fi, the map which connects the velocities from just be- 
fore the kick to the next kick (by the use of Eqs. (|3l|)- 



^ = 1- 



n + l 



T 

7+ to 



dx\ O N fdy\ , n , 
COS—T + -| sin( T) 
dtj + 7+ \dtj + 7+ 



T 



(38) 



rt+l 



7+ TO 



dx\ fd y \ , n . 

~ sin( — T) + cos — T) 

dt J, S+ ' \dt J. S+ ' 



(39) 



, _7+ 

n 

7+ 
O 2 r 



n.+ l 



dx\ 

~dt) n+1 



dy 
dt 

dx 
~dt 



exact solution 



semi-exact 
solution 



approximated 
solution 




numerical 
integration 

by elementary 
Fresnel functions 

exact 
solution 



(40) 



Hi 



Notice that the minus signs in Eq. (37) refer to the n 
points of the map. 

It is possible to return to the non-radiative limit by 
letting i — * 0. In that case the map is equivalent to the 
map which was derived by Longcope and Sudan || . The 
nonrclativistic limit, which was derived by Zaslavskii et 
al jjj, is achieved by letting ^- — > and c — > oo. 



FIG. 2. A diagram which present the all possibilties for 
the map construction. The wider line indicates the chosen 
combination. 



The radiative map is then 



(dx\ 
\ dt)- 



(*). 



sin a 



n+i cosh — + \/l — 81 sinh — 

T V + T ° 



(41) 
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'dy\ 

Xn+l 



(f ) + sina 



(*). 



cos a 



cosh — 

TO 



It 

7+ 
Wt 



1-01 sinh -t 



n+i v * / + 
/ dx 



(42) 



(43) 



where the values immediately after the kick are as in Eqs. 
© and a is defined in Eq. ( p7| ) (a — > —a since e = — eo) 
. Note that for our limit {T/tq -C 1) there is not a great 
difference between the map given by Eqs. (pq)-(|40[) and 
Eqs. ©-dD. 



V. RESULTS 

In order to obtain a web structure, there are two con- 
ditions that have to fulfilled. In the nonrelativistic case, 
the first condition is that the ratio between the gyration 
frequency, f2, and the kicking time, T, is a rational num- 
ber, a condition which can be expressed as follows, 



P 



(44) 



where p, and q, are integer numbers. Secondly, one must 
start from the neighborhood of the unstable fixed point 
(otherwise, the particle will not diffuse, and will not cre- 
ate a web structure), 



=r) =0 



(2n+l) — , 



(45) 



where n is an integer number. The symmetry of the 
web is determined by p and q. If, for example p = 1 
and q = 4, the particle is kicked four times during one 
cycle, and thus, the symmetry of the web will be a four 
symmetry ffl. 

In the relativistic case [0, the above conditions 
are slightly different, because of the additional factor 
\Jl — 0q which multiplies the OT term. However, if the 
initial velocities are small, i.e. vo <C c, the additional 
factor is close to 1, and thus the structure of the web 
should not change. In the radiative and relativistic case, 
as well as the non-radiative relativistic case, conditions 
(Q) and (f4|) become, 



Pi =o 

dtJo 
' dy 

dt 



(2n+l)Jl-0l 



(46) 



For sufficiently large initial velocities the above condi- 
tions do not hold anymore, and the web structure is not 
observed. 

In this section we will compare (qualitatively; a quan- 
titative treatment for the diffusions rate will be studied 
elsewhere) the diffusion of the non-radiative particle and 
the radiative particle. Intuitively, one would expect that 
a radiative particle will diffuse more slowly than a non- 
radiative one, since the radiation effects act like friction, 
and are thus expected to "stop" the particle. However, 
this naive explanation is not true, as we will demonstrate 
below. 

In order to investigate the above assumption, we have 
chosen a four symmetry structure (q = 4). We used 
the same initial conditions for all cases; just the kick- 
ing strength K has been changed. The initial conditions 
were B — 10, v x $ — 0, v y fl = 4 x 10~ 5 , k = 10, and 
the number of iterations was N — I0 6 . In Fig. ^, we 
present the results as follows: the first column is the non- 
radiative case, the second column is the radiative case, 
and the third column is the total velocity versus the it- 
eration number (the circles indicates the non-radiative 
case and the squares the radiative one). In the first row 
K = 1 x I0~ 5 , in the second K = 2 x 10~ 5 , and in the 
third K = 4 x 10~ 5 . 

As can be seen clearly from Fig. ||, the web structure is 
valid for the radiative case (top panel, first column). The 
web structure is also exists in the other cases; the web 
cells are very small and they are difficult to see. Although 
the radiative particle and the non-radiative particle start 
from the same initial conditions they behave differently. 

In some cases the diffusion rate of the radiative case is 
larger then the diffusion rate of the non-radiative case, for 
example, in the top and the bottom panels after 9 x 10 5 
iterations the radiative particle reaches a higher velocity 
then the non-radiative one. On the other hand, in the 
middle panel the radiative particle is slower then non- 
radiative one. Obviously, it is impossible to draw any 
general conclusion for the question of whether a non- 
radiative particle is faster then the radiative one. For 
this a more systematic approach should be carried out. 
This would be beyond the scope of the present paper and 
will be considered elsewhere. 

When the velocity of the particle reaches close to the 
velocity of light it actually stops growing (according to 
evolution in the time t) although it increases its energy. 
In Fig. ^ we show an example for that behavior. The pa- 
rameters values are : B = 10, v X: o = 0, v y fl = 1 x 10~ 5 , 
k = 10, K = 0.01, and the number of iterations was 
N = 1 X 10 6 iterations (every 10 th iteration was plotted). 
Fig. Ha and Fig. Ep are the non-radiative cases, while 
Fig. |p and Fig. p are the radiative cases; in Fig. ^a 
and Fig. the map of v x versus v y is presented, while in 
Fig. |l|c and Fig. ^d the total velocity v versus the itera- 
tion number is plotted (every 10 th iteration was plotted). 
As seen from Fig. EL the particle accelerates quickly to 



G 



high velocity, and most of the time stays with this high 
velocity, approaching light velocity asymptotically, while 
increasing its energy to infinity. 

VI. SUMMARY AND DISCUSSION 

In the present paper we have investigated the effect of 
radiation on the stochastic web. Under the restriction of 
small magnetic fields (B < O(IOOT)) an iterative map 
was constructed. Moreover, the effect of radiation on the 
stochastic web is very small because of the small magnetic 
field. Qualitatively the non-radiative and the radiative 
cases have similar web structure. Despite of the naive 
expectation that the diffusion rate in the radiative case 
should be smaller then non-radiative case, one can find 
cases in which the opposite effect is observed. 

Although it seems the effect of radiation is not quali- 
tatively significant under laboratory conditions, it might 
have a strong effect in the presence of a strong magnetic 
field. Such a magnetic field could occur near or in a neu- 
ron star (or other heavy stars) and it can cause a large 
radiation correction to the motion of the particle. In that 
case one does not expect to have a web structure even 
when the particle has very small velocity, since the par- 
ticle might move from the unstable fixed point and thus 
destroy the web structure. On the other hand, another 
scenario can take place; the magnetic field may be so 
strong that the particle almost stops its motion before 
the next kick, and under some symmetry conditions it 
can diffuse to the large energy regime. The described 
scheme can reflect new phenomenon not previously pre- 
dicted, and could give some insight to the stellar systems. 

In the present study we did not consider the radia- 
tion during the S kick because of lack of proper technical 
methods as well as analytical knowledge. Taking into 
account the radiation during the kick could change our 
conclusions, and as for other open questions which men- 
tioned in this paper, this point will be treated elsewhere. 
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Fig. 3 
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FIG. 3. Different types of diffusion behavior as described 
in the context. 



Fig. 4 




FIG. 4. The behavior of a particle with high velocity, (a) 
the nonradiative map map, (b) the radiative map, (c) the total 
velocity of the non-radiative case, and (d) the total velocity 
of the radiative case. 
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